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The time evolution of thermal states of a mirror released from a tight harmonic trap is studied. 
After the release no dissipation is assumed to be present and the mirror is, after a time delay, 
kicked into a momentum superposition state. The thermal character of the initial state washes out 
the telltale interference patterns of the superposition but no loss of coherence is found. This in- 
vestigation resolves a controversy about decoherence-without-dissipation and shows that entrained 
measurements can be surprisingly insensitive to temperature effects. 
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I. INTRODUCTION 

For massive objects it is hard to create a pure quantum 
state unentangled with the environment. Many studies 
of the decohering effects of the environment do however 
assume that the studied system is in a pure state initially, 
unentangled with the environment. Although realizable 
for atomic systems this is implausible for massive objects 
such as opto-mechanical mirrors pQ. Ford et al. 00] are 
therefore quite justified in insisting that for investigations 
of the dynamics of quantum superposition states of mas- 
sive systems, weakly coupled to the environment, thermal 
initial states should be considered. They derived expres- 
sions for the loss of interference contrast of formally ther- 
malized momentum superpostion states of free massive 
objects. Their interpretations have received a somewhat 
mixed reception though. Notably, their interpretation 
of a process they dubbed "decoherence without dissipa- 
tion" [3J has been either criticised E] or ignored [5] . 

Here, following Ford et al. 's approach, a mirror is stud- 
ied which initially is weakly coupled to the environment 
such that we can neglect dissipation other than assuming 
the mirror is in thermal equilibrium with the laser beams 
that form a tight harmonic potential holding and cool- 
ing it [9 . When the lasers are switched off the mirror 
becomes effectively decoupled from the environment, we 
model this as a ballistically expanding mirror without 
dissipation 00] ■ Once the coherent width of the mir- 
ror's spatial distribution is sufficient to imprint a pho- 
ton's momentum-recoil pattern the mirror is kicked into 
a momentum superposition state and then investigated 
using the 'entrainment procedure' of reference [10] . 

This work addresses two main questions. Can deco- 
herence without dissipation occur? How badly does 
the thermal nature of the initial state affect the forma- 
tion and persistence of interference patterns of a mirror's 
wave packet [TU]? 

Ford and O'Connell use the attenuation of the inter- 
ference contrast as a measure for decoherence. This is 
not the conventional approach for the quantification of 
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decoherence, which is based on the loss of off-diagonal 
elements of the reduced density matrix .5 8 . One can 
therefore refute Ford and O'Connell's assertion, that they 
found evidence for "decoherence without dissipation" , 
on formal grounds, arguing, that no such process can oc- 
cur without coupling to the environment [3 IS] • I tem- 
porarily adopt Ford and O'Connell's approach though 
and investigate the visibility of the interference pattern, 
since it is the quantity an experimentalist is most likely 
to measure . 

Ford et al. 00] consider a free particle in a Gaus- 
sian state with arbitrary width which is then formally 
thermalized by convolution with a Maxwell-Boltzmann 
velocity distribution function. One could argue, using 
Ford, Lewis, and O'Connell's own logic 0, that the ar- 
bitrary spatial extension and precise form of the initial 
Gaussian state on the one hand, and its formal ther- 
malization at an arbitrarily chosen temperature 6 on the 
other, might not be independent and the details need 
further justification. Indeed, the environmental tempera- 
ture typically determines the spatial extent of emerging 
position-pointer states and a generic mechanism, such 
as photons randomly scattering off point particles with 
spatially extended center-of-mass wave function, can lead 
to wave functions that do not have Gaussian shape |llj . 



To avoid such ambiguities a possible physical imple- 
mentation of the scenario modelled here is introduced 
in section |TT] Section Mil illustrates different cases of our 
scenario employing varied physical parameters. In sec- 
tion |IV| we consider the impact of temperature on the 
expected visibility of the interference pattern and its loss 
(and resurgence) over time. Finally, in section[VJ we con- 
sider the 'free case', mostly discussed by Ford et al., we 
will see that it can be mapped onto the case of a mir- 
ror trapped harmonically. A harmonically trapped ther- 
mal mirror shows periodic resurgence of the interference 
pattern, we therefore conclude that "decoherence with- 
out dissipation" is absent even when adopting Ford and 
O'Connell's approach using interference contrasts instead 
of non-diagonal density matrix elements. 
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II. STEPS OF THE SCENARIO 



state' (it is larger for 'released excited states' 



A. Trapped and cooled mirror 

A mirror with mass M, initially trapped by a stiff har- 
monic potential with spring constant Kq and cooled to 
temperature 9 [5], is described by the thermal density 
matrix 
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For M = 10 
stants Kq » 
60 [im-s^ 1 to 
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to 10 
10 6 Nm" 
10 nm-s 



~ 10 kg pQ and optical spring con- 
_1 [9] this velocity ranges from 
. In other words, for a mirror 



mass of 10 _15 kg, the free ground state wavefunction ex- 
pands to an optical wavelength size of roughly 500 nm 
within an expansion time of approximately 9 ms. 



Here 6^ = hy/K /M jk B is the associated Einstein- 
temperature where h is Planck's and ks Boltzmann's 
constant. In all numerical calculations we will set H = 1 
and ks = 1- For opto- mechanical oscillator masses rang- 
ing from M = 10~ 15 to 10~ 10 kg [1J and recently reported 
large optical spring constants Ko s» 10 6 Nm _1 [5] the 
oscillator's resonance frequency Sl = \J Kq/M ranges 
from 3 • 10 10 s _1 to 10 8 s _1 and the associated Einstein- 
temperature Qe = Mlo/kB from 0.24 K to 0.76 mK. 
Such temperatures can be reached in optical cooling se- 
tups HH H3] . We will see that the scheme presented 
here can tolerate temperatures up to several Qe of the 
trapping potential. 

The normalized harmonic oscillator eigenfunctions 
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meterized by the system's ground state position standard 
deviation 



(2) 



For the values for Kq and M cited above, this position 
spread is very small {ctq ps 10~ 14 to 5 • 10~ 16 m). This 
is good to protect the oscillator from decoherence but 
makes it hard to see interference fringes when momentum 
superposition states are formed [10] . 



B. Released mirror 

To be able to 'imprint' the interference pattern asso- 
ciated with the momentum superposition state |10j we 
will therefore assume that the stiff trapping potential is 
suddenly switched off so that the mirror is very weakly 
trapped or set free. In the free case, the ensuing ballistic 
expansion of the stiff potential's groundstate yields the 
'free groundstate' wave function [14] 
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is the time-dependent position spread which gives rise to 
the ballistic expansion velocity of the 'released ground- 



C. Kicked thermal mirror 

In order to achieve the momentum superposition, it is 
suggested to use an interferometric setup featuring a half- 
silvered beam splitter that splits a single photon (gen- 
erated in spontaneous parametric down-conversion |15j ) 
into two equally strong partial waves. These partial 
waves are directed onto either side of the mirror and, 
upon reflection, transfer their momentum recoil to the 
mirror in a desirable fashion: namely, after tracing over 
the photon wave function the mirror's wave function be- 
comes multiplied with an effective kick-factor |10j 



K, = — i sm(2nx — — ) , 
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where the relative phase <f> is fixed by the interferometric 
setup and the port the photon takes [10]. Here, it is 
assumed that the mirror has unit reflectivity and k is 
the photon's effective wavenumber, associated with its 
normal momentum component, for more details and the 
case of imperfect mirrors see [10] . 

The 'released eigenstates', spawned by the 
eigenfunctions ip n of the stiff potential, are la- 
belled W n (x, t;0,T, k, Ko); here r < denotes the 
time when the mirror is released from the stiff potential 
K x 2 /2 into a weak potential kx 2 /2 (or set free: k = 0, 
in which case * (^, 0; 0, -t, 0, K ) = $ (x, t) of Eq. 

Then, at time t = 0, the photon kicks the mirror into 
a momentum-superposition state. Through the momen- 
tum kick the wave functions become combined to form 
the desired momentum superposition states T ra oc IOS n 

T n (x, t;p 7 , t, <fi, k, K ) = J\f n [^ n (x, t;p^,r, k, Kq) 

- e^ n {x,t:~ Pl ,T,k,Kv)]. (7) 

Here, M n is a normalization constant that absorbs over- 
all phase factors and primarily depends on the size 
of the photons' momentum transfers p 7 = 2hn, com- 
pare Eq. ([6]), and the position spread of each state ^„ at 
the time (t=0) when the photon-kick occurs. The ther- 
mal momentum superposition state thus has the form 
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FIG. 1. Probability densities P(x,t) of quantum harmonic 
oscillator with mass M = 1 released into a weak harmonic 
potential (r = —T/4, k — 1) and subsequently momentum- 
kicked at time zero (<f> = n, for how fixed phases arise in 
entrained measurements see QH]). Plots A and B show re- 
leased and kicked (p 7 = 1) groundstate with different 
initial conditions [Kq = 2 4 and single kick, versus, Kq = 3 4 
and momentum superposition kick). Similarly, plots C and 
D show released and kicked (p 7 = 2.5) third excited state Vf^ 
with K = 2 4 and r = -T/4, and K = 3 4 and r = -T/8, 
respectively. 

The states T„ up to 13-th order have been determined in 
this work, so, the cut-off N — 13. Upon inclusion of the 
'free' kicked groundstate T , these lowest fourteen oscil- 
lator states allow us to model the stiff oscillator at tem- 
peratures of up to three times the Einstein-temperature: 
Er=i4 e " /3 ° < °- 01 ' since this accounts for 99% of the 
probability distribution for temperature 6 — 3.0 • Qe at 
which a harmonic oscillator has roughly 99.08% of the 
classical heat capacity |fcs. At such temperatures, the 
character of the system can be considered to be between 
quantum and classical. 



III. DYNAMICS OF PURE STATES 

We now study the dynamics of the groundstate, which 
is an even function in position, and the third excited 
state, which is odd. 

Figure [I] A shows the groundstate released from the 
stiff oscillator potential at r = — 3T/4 where T — 
2TT V / M/k = 2tt/ u> is the mirror's period time. It plays 
out as an underlying squeezed and anti-squeezed breath- 
ing motion in evidence in all plots of Figure [T] At time 
zero the mirror is kicked to the right and so shows a com- 
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FIG. 2. Probability densities P(x,t) = p(x,t; §,-§,0,ifo) of 
quantum harmonic oscillator with mass M = 1 set free (k = 
0) at time r = —1/2 from stiff oscillator with spring con- 
stant Ko = 2 4 and momentum-superposition-kicked at time 
zero: p 7 = 3/5. Plot A shows mirror cooled almost to ground 
state 9 = 0b/3. It demonstrates that the small extension of 
the wave packet is not wide enough to imprint a full interfer- 
ence pattern (</> = ir), although a slimming and heightening 
of the probability distribution can be seen. Plot B same as A 
for 9 = 3Gb. Clearly the wave packet is wider and expands 
faster at this higher temperature allowing for a full imprint 
of the interference pattern. Plots C, same as B with cen- 
tered interference minimum (</> — 0). The dephasing of the 
interference pattern in this case is greater than in B because 
fewer states contribute to dephasing on the fringe than at the 
center x = 0. Plot D: Under the dotted envelope P(x, 0) the 
thermally weighted contributions of the 'released eigenstates' 
are shown, demonstrating that fewer eigenstates contribute to 
interference on the fringe (case B) than at the center (case C). 



bination of breathing and oscillation. Plot B shows the 
same scenario except for the facts that the release from 
the stiff potential happens at time r = —T/4 and the 
momentum-kick is in superposition. The formation of 
the superposition is highlighted by the red line immedi- 
ately before the momentum kick and the yellow line when 
it happens. The same yellow profile line has been moved 
downstream and shows that intensity patterns re-form 
every half period. 

Plots [T] C and D show analogous scenarios for the 
third excited state released from the stiff potential at 
times t = —T/4 and —T/8 respectively. Plots B and D 
illustrate that the interference pattern of the momentum 
superposition is imprinted onto the underlying spatial 
probability distribution of the state with full contrast. 
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IV. DYNAMICS OF THERMAL STATES 



A. Benchmark for dephasing function 



At the temperatures considered here, ranging from 
zero up to three times the Einstein temperature, the in- 
terference with full contrast can be observed for all set- 
tings of the system's parameters such as spring constants, 
mass, and photon momenta, provided that the released 
mirror wave functions had time to ballistically expand 
sufficiently widely to accommodate an imprint at the ef- 
fective wavelength [TO] of the photon recoil when it is be- 
ing kicked into the superposition state, compare Fig. [2] A 
versus B. In other words, wave packet expansion and the 
momentum-superposition interference imprinting opera- 
tions do not commute. Each wave function T n is by itself 
fully coherent and the momentum superposition imprints 
occur at the same spatial positions for all of them. So, 
the effects add up to an interference pattern with full con- 
trast. This implies that incoherent mixtures consisting of 
states spread out wide enough to carry the momentum 
superposition imprint can be endowed with an interfer- 
ence pattern when using the entrainment procedure of 
reference [TU]. This observation, which also applies to 
states with temperatures much higher than studied here, 
is one of the central messages of this work. 

Note that the dephasing of the imprinted interference 
pattern does not only depend on the temperature but also 
its location. Further out from the center fewer 'released 
eigenstates' T„ contribute to the thermal wave packet 
and therefore dephasing is slowed down there as com- 
pared to dephasing at the center, x = 0, compare Fig.[2]B 
with C. On the fringes the interference gets washed out 
more slowly. 

Additionally, more highly excited 'released eigenstates' 
have a wider spread and higher ballistic expansion ve- 
locity, they therefore can carry an interference imprint, 
when lower lying 'released eigenstates' might not be able 
to do this, see Fig. [2] B as opposed to A. Surprisingly, 
higher temperatures can help in generating an interfer- 
ence pattern imprint because they expand the width of 
the mirror's wavepacket. 

Since they have different energies the various wave 
functions T„ evolve at different rates and therefore mu- 
tually dephase thus washing out the interference imprint. 
The general scenario is varied and its detailed quantifica- 
tion is beyond the level of this work. We would, however, 
like to point out that the breathing motion underlying the 
dynamics in the trapped case (see breathing motion of 
pure states in Fig. [T]) is least pronounced at every turn- 
ing point or odd quarter period point: (2n + l)T/4 = 
(2n + l)27r/(4w) counted from the moment of release r. 
This is confirmed by Fig.[3]B and C which show slowed 
dephasing since the momentum-superposition kick hap- 
pens at the quarter time (r = —T/4, t = 0). 



The dephasing times of the interference patterns in evi- 
dence in Figs. [2] and [3] can be crudely benchmarked. To 
do this we derive a reference approximation A. In the 
thermal state the wave functions are incoherently added 
up and so the dephasing is due to the summing up of 
their individual dynamics, weighted with the Boltzmann- 
factors. To derive our benchmark approximation we form 
a coherent instead of an incoherent sum of wave func- 
tions using the dominant differences in their energy ex- 
pressions, namely the harmonic oscillator energies of the 
initial stiff trapping potential, in the time evolution phase 
factors. This neglects variations due to the transferred 
kinetic energies but should be compensated for by the 
overestimate due to the formation of a coherent sum. 
Moreover, we assume that all wave functions contribute 
equally (weight '1') rather than with their true local 
weight |T„| 2 , this also overestimates the dephasing effect. 
This yields an expression containing the 'Boltzmann- 



sums' b(0,t) = £~ =0 e 
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Since the probability distribution of the released mir- 
ror does not oscillate at frequency but at twice the 
frequency w of the weak trap, it is released into, our ap- 
proximation for the interference pattern visibility adopts 
this and reads 



A(0,t) 
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Clearly this approximation is inconsistent; for instance 
the probabilities instead of the amplitudes of the 
Boltzmann-factors are used, but then again A(9/2, t) 
does use the amplitudes, in short: the approximation 
A can only serve as a guide to the eye, compare Fig. |4j 
but it allows us to underline the fact that dephasing hap- 
pens on time scales of the system's evolution ~ l/(2w). 
This answers the second question posed in the introduc- 
tion: Since the approach using entrained measurements 
allows for very fast state preparation and analysis [TU] , we 
can hope that, assuming it can be experimentally imple- 
mented, thermal dephasing does not affect it too badly. 



V. NO DECOHERENCE WITHOUT 
DISSIPATION 

The state ^ shows full interference contrast at at the 
time (t = 0) when it is kicked into the momentum su- 
perposition state [TU] . The associated visibility V of the 
mirror's interference pattern at the origin, as a function 
of time, is given by the expression 

_ \p(Q, t;p 7 , t, 0, k, K ) - p(0, t;pj, r, tt, k, K )\ 
p(Q, t; p 1 , r, 0, k, K ) + p(0, t; p 1 , r, tt, k, K ) 



FIG. 3. Probability densities P(x,t) of quantum harmonic 
oscillator with mass M = 1 released from stiff oscillator po- 
tential into weak potential and kicked into momentum super- 
position state at time t — 0: Plot A 6 = 30e, K — 16, 
k = 1, p 7 = 2 and r = -T/20; B = 36b, # = 64, fe = 4, 
p 7 = 2 and r = -T/4; C 9 = 2.8159b, if = 32, fc = 2, 
p 7 = 1 and r = -T/10; D N 38s, if = 16, k = 1, p 7 = 4 
and r = —T/4. It can be seen that the interference patterns 
always re-form with full contrast after half a period. More- 
over, the time over which the interference pattern persisters 
is lengthened when one kicks the mirror at the point of great- 
est expansion rather than earlier (see B as opposed to A). 
The dephasing of the interference pattern is weaker for longer 
wavelength imprints (smaller momenta), compare A and B, 
with C and D. 



It has this definite form -without the need to maximize 
or minimize over the settings of the phase </>- because 
all constituent wave functions have definite parity in the 
^-coordinate. 

The visibility degrades over time because of dephas- 
ing of the components of the thermal density matrix (|8| . 
But, if the mirror is released into a weak harmonic trap- 
ping potential with spring constant k instead of being 
set free (fc = 0), its wave function will re-form the same 
probability distribution (or its position-inverted mirror 
image) every half cycle of length T/2; this is illustrated 
in Figs. [T] and [3] 

A. Map between free and trapped states 

In theory, the ballistically expanding wave functions 
can be determined explicitly, using suitable propagators. 
In practice, this turns out to be surprisingly difficult 
even when using a contemporary computer algebra sys- 



FIG. 4. Same parameters as Fig. [3] (except for temperature 9): 
plots of the decadic logarithm ld\V) of the visibility V as de- 
fined in Eq. [TT] as a function of temperature and time (col- 
ored mesh) together with the benchmark approximations A 
of Eq. 10 (brown contour-sheets). This figure shows how with 



increasing temperature the visibility of the interference pat- 
tern washes out faster but always reforms fully after half a 
cycle. Moreover, if the mirror gets kicked into the momen- 
tum superposition state near the turning point (r = —T/4) 
the interference pattern persists longest, compare text and 
Fig. |H 



tern [16]. An eigenstate of the initial stiff potential has 
overlap with many states of a weak trapping potential 
it is released into, making the resulting expressions un- 
wieldy. Instead, a trick is employed here. A mapping has 
been used to map between harmonically trapped and free 
states, for a detailed discussion see [T3]. This way it be- 
came possible to generate Figs. [3] and [4] 

Aside from being a mathematical trick this mapping 
also shows that no decoherence occurs for the case of the 
free particle, as was claimed in references [2R3]- Instead, 
the wave functions dephase but the mapping back to the 
trapped case and further propagation in time shows that 
the interference pattern reforms with full contrast. Note 
that, although challenging in practice for a widely spread 
out wave packet, in principle the mapping corresponds 
to the experimentalist switching on the trapping poten- 
tial [14]. This answers the first question posed in the 
introduction: there is no decoherence without dissipa- 
tion. 

It raises the new question as to when to let the pho- 
ton impinge on the mirror to imprint the superposition 
pattern. The answer to this question depends on the 
goals one wants to achieve experimentally. If one waits 
for time T/2 after the release the wave function has 
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expanded furthest. This stage maybe be advantageous 
if the large mass of the mirror makes it impossible to 
achieve sufficient position spread early in its evolution. 
On the other hand, imprinting the interference pattern 
at this stage implies that the momentum kick is trans- 
ferred to the mirror at a time of largest spatial exten- 
sion and subsequently, the mirror will contract again, see 
Fig. [T]B. Unless very large recoil values are a employed 
this precludes the formation of separate spatial humps 
in the probability distribution which are typically con- 
sidered to be the hallmark of so-called "Schrbdinger-cat" 
spatial superposition states. 



a very weak potential (or set free) [T7J HH] . It has shown 
that decoherence without dissipation, based on thermal- 
ized initial states alone, does not occur as envisaged in 
references 2 -4 . It also studied the formation and persis- 
tence of interference patterns of the mirror's wave packet 
after it is kicked into a superposition state. It was shown 
that dephasing and re-formation of the interference pat- 
terns happens on the time scale of the system's time evo- 
lution and should therefore be detectable when using suf- 
ficiently fast detection schemes such as the 'entrainment 
procedure' sketched in reference [TU] , 



VI. CONCLUSION 



This work studied effects of thermalized initial states 
of a harmonically trapped quantum mirror released into 
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